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Abstract 

A key problem in multiobjective linear programming is to find the set of all efficient 
extreme points in objective space. In this paper we introduce oriented projective geometry 
as an efficient and effective framework for solving this problem. The key advantage of 
oriented projective geometry is that we can work with an "optimally simple" but unbounded 
efficiency-equivalent polyhedron, yet apply to it the familiar theory and algorithms that are 
traditionally restricted to bounded polytopes. We apply these techniques to Benson's outer 
approximation algorithm, using oriented projective geometry to remove an exponentially 
large complexity from the algorithm and thereby remove a significant burden from the 
running time. 
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1 Introduction 

In traditional linear programming, our goal is to maximise a linear objective function over some 
polytope X C E„ (where E„ denotes n-dimensional Euclidean space). In multiobjective linear 
programming, we attempt to maximise many linear objective functions over X simultaneously. 
Multiobjective linear programming is an important practical tool in multicriteria decision mak- 
ing [ini [Ml [33 [S], and recent decades have seen a surge of activity in understanding the 
complexities of such problems and developing algorithmic tools to solve them. 

Typically different objectives achieve their global maxima at different points of the polytope 
X, resulting in trade-off situations in which some objectives can be improved at the expense 
of others. In the most general solution to a multiobjective linear program, we seek to find all 
efficient (or Pareto- optimal) points, which are those solutions from which we cannot increase 
any one objective without simultaneously decreasing another (21j . 

The set of all efficient points in X is difficult to manage, since it is typically large and non- 
convex. Instead it is preferable to seek a smaller combinatorial description of the efficient set. 
A common approach (which we follow in this paper) is to compute the efficient extreme points 
of X (that is, the efficient vertices of X), from which the full efficient set can be reconstructed 
[T1[I1[15]. 

Early approaches to multiobjective linear programming were based on techniques from the 
simplex method [3 [TTJ [301 [S]; and involved walking through the polytope X. However, for 



problems with large numbers of variables n, the polytope X C E„ can grow to become extremely 
complex. 

A new approach was therefore developed to work in the smaller objective space Ep, where 
p is the number of objectives under consideration |121 [Ml [15] . The original polytope X C E„ 
naturally maps to an outcome set C Ep, which is a polytope that describes all achievable 
combinations of objectives. 

Working in objective space is preferable because the running time of many relevant algorithms 
is exponential in the dimension of the polytope |36[|42) — this means that working in p dimensions 
instead of n can have an enormous impact on the computational resources required. Moreover, 
it is often easier for the decision maker to choose a preferred outcome from instead of a 
specific combination of input variables from X [SI I12j . 

We therefore focus on the problem of finding all efficient extreme points of the outcome set 
Y^. There are many algorithms in the literature for solving this problem. One notable family 
of algorithms is based on weight space decompositions [H [HI [30] > extending earlier theoretical 
work on the vector maximum problem [24i 140] . Another important family of algorithms is 
based on outer approximation techniques from the global optimisation literature [3J [5] . It 
remains relatively unknown how these algorithms compare in practice, and it is possible that 
both families have important roles to play as the problem size becomes large [Hj. 

A key optimisation in these algorithms comes from working with ejjiciency-equivalent poly- 
hedra. These are polytopes and polyhedra in Ep whose efficient points are the same as for the 
outcome set Y^ , but whose non-efficient regions are much simpler. Two such polyhedra feature 
prominently in the literature: 

• Dauer and Saleh [12|, 113] define the efficiency-equivalent polyhedron Y-, which is combi- 
natorially simple but which has the practical disadvantage of being unbounded (extending 
infinitely far in the negative directions). Being unbounded means that it cannot be ex- 
pressed as the convex hull of its vertices, which causes problems for important procedures 
such as the double description method of Motzkin et al. [29] . 

• Benson [21 [3l [5] describes the efficiency-equivalent polytope Y^, which is bounded with 
a cube-like structure. Although it is less combinatorially simple, the bounded nature of 
Y^ makes it more practical for use in algorithms, and indeed Y^ plays a key role in the 
family of outer approximation algorithms mentioned above. 

In this paper we develop a new framework for algorithms in multiobjective linear program- 
ming, based on oriented projective geometry. Introduced by Stolfi [37], oriented projective ge- 
ometry augments Euclidean space with additional points at infinity. Unlike classical projective 
geometry, oriented projective geometry also preserves notions such as convexity, half-spaces and 
polytopes, all of which are crucial for linear programming. 

We show that by working in oriented projective p-space, we can achieve the best of both 
worlds with efficiency-equivalent polyhedra. Specifically, we work with an efficiency-equivalent 
projective polytope T-, which inherits both the simple combinatorial structure of Y- and the 
wider practical usability of y^. 

To show the efficiency of this new framework, we undertake a combinatorial analysis of 
the bounded Euclidean polytope Y^ and compare this to the projective polytope T-. We 
show that the non-efficient structure (that is, the "uninteresting" region) of Y^ is exponentially 
complex: the number of non-efficient vertices is always at least 2^ — 1, and can grow as fast 
as ri(nL(P~i)/^J ). In contrast, we show that the projective polytope T- has just p non-efficient 
vertices, a figure that is not only small but also independent of the problem size n. 
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Even for problems where the number of objectives p is smah, our framework enjoys a sig- 
nificant advantage by using a projective polytope whose non-efficient region has a complexity 
independent of n. For problems where p grows larger, such as the aircraft control problem of 
Schy and Giesy with 70 objectives [32], or applications in welfare economics where p can grow 
arbitrarily large [3T] , even the lower bound of 2^ — 1 puts at a severe disadvantage compared 
to the projective polytope T- that we use in our framework. 

To show the effectiveness of this framework, we apply it to the outer approximation algorithm 
of Benson [5]. The result is a new algorithm that works in oriented projective space, combining 
the practical advantages of the original algorithm with the combinatorial simplicity of T-. 
By studying the inner workings of this algorithm, we show how the complexity results above 
translate into significant benefits in terms of running time. 

The layout of this paper is as follows. Section [2] summarises relevant definitions and results 
from polytope theory and multiobjective linear programming, and Section [3] offers a detailed 
introduction to oriented projective geometry. In Section 0] we discuss efficiency-equivalent poly- 
hedra, and prove the combinatorial complexity results outlined above. We apply our techniques 
to the outer approximation algorithm in Section[Sl and in SectionlHlwe conclude with a discussion 
of the theoretical and practical implications of our work. 

2 Preliminaries 

In this brief section, we outline key definitions from the theory of polytopes and polyhedra and 
the theory of multiobjective linear programming. For polytopes and polyhedra we follow the ter- 
minology of Ziegler [32] and Matousek and Gartner [35]. For multiobjective linear programming 
we follow the conventions and notation used by Benson [5]. 

Throughout this paper we work in both traditional Euclidean geometry and the alternative 
oriented projective geometry. We denote d-dimcnsional Euclidean space by E^, and each point 
in Erf is described in the usual way by a coordinate vector x — [xi, . . . , Xd) E K'*. Although the 
geometry E^ and the coordinate space R"* are identical in the Euclidean case, this distinction 
becomes important in oriented projective geometry, where we describe c?-dimensional geometry 
using coordinates in We delay any further discussion of oriented projective geometry until 

Section [3] 

The concepts of polytopes and polyhedra are fundamental to linear programming. A poly- 
hedron in Ed is an intersection of finitely many closed half-spaces in E^; this may be either 
bounded (like a cube) or unbounded (like an infinite cone). A polyhedron that is bounded is 
also called a polytope. Every polytope in can be expressed as a convex hull of finitely many 
vertices in E^, whereas unbounded polyhedra cannot be expressed in this way. 

A hyperplane in E^ is a (d — l)-dimensional affine subspace, and consists of all points with 
coordinates {x G K"^ | h • x = a} for some h € R'^ and a G R. If P C E^ is a polyhedron 
(bounded or unbounded) and iif C E^ is a hyperplane, then we call H a supporting hyperplane 
for P if (i) H contains some point of P, and (ii) all of the points in P lie within H and/or to 
the same side of H . In other words, (i) h • x = a for some x G P, and (ii) either h • x > a for 
all x G P, or else h • x < a for all x G P. The intersection of P with a supporting hyperplane 
is called a face of P0 

^For theoretical convenience, two additional faces are defined: a "full face" containing all of P, and an "empty 
face" containing no points at all. However, neither of these is relevant to this paper. 
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The dimension of a polyhedron is defined to be the dimension of its affine hull (in particular, 
a polyhedron in must have dimension < d). Likewise, the dimension of a face is defined to 
be the dimension of its affine hull. If P is a polyhedron of dimension k, then faces of dimension 
0, 1 and {k — 1) are called vertices, edges and facets of P respectively. 

A multiobjective linear program requires us to maximise p linear objective functions over 
a polyhedron in E„. This polyhedron is called the feasible solution set X, and consists of all 
X £ E„ for which Ax = b and x > 0, where A is an m x n matrix and b e M™. Like Benson 
[S], we assume for convenience that X is bounded; that is, X is a polytope. 

The p linear objectives are described by a p x n objective matrix C, and our task is to 
maximise Cx over all x G X. We define the outcome set as {Cx | x e X}. The outcome set 
is in turn a polytope in Ep, and our task can be restated as maximising y e Y^. 

In general we cannot achieve a global maximum for all p objectives simultaneously, so instead 
we focus on non-dominated outcomes. We call an outcome y € Y^ efficient (also called Pareto- 
optimal or non- dominated) if there is no other y' G Y^ for which y' > y. The efficient outcome 
set Y^ is defined to be the set of all efficient outcomes in Y^ . 

More generally, an efficient point in any polytope or polyhedron Z C Ep is a point z G Z for 
which there is no other z' G Z satisfying z' > z. The set of all efficient points in Z is denoted 
Ze, and this set can always be expressed as a union of faces of Z. If Ze = Y^ then we call Z 
an efficiency- equivalent polyhedron for Y^ . 

Since the efficient outcome set Y^ is generally infinite, we focus our attention on the set 
of efficient extreme outcomes. An extreme point of a polyhedron is simply a vertex of the 
polyhedron, and an efficient extreme outcome is an efficient vertex of Y^ . Because Y^ is a 
polyhedron there are only finitely many efficient extreme outcomes, and from these we can 
generate all of Ye if we so desire [5D]. Our focus then for the remainder of this paper is on 
generating the finite set of efficient extreme outcomes. 

3 Oriented projective geometry 

In classical projective geometry, we augment the Euclidean space E^ with additional "points at 
infinity" . This offers significant advantages over traditional Euclidean geometry. Theoretically, 
we achieve a powerful duality between points and hyperplanes. Computationally, we can perform 
simple arithmetic on infinite limits (such as the "endpoints" of lines and rays); moreover, we 
lose the notion of "parallelness" and thereby avoid a myriad of special cases. For a thorough 
overview of classical projective geometry, the reader is referred to texts such as Coxeter or 
Beutelspacher and Rosenbaum [7^. 

A key drawback of classical projective geometry is that we also lose important concepts 
such as convexity, half-spaces and polytopes. Oriented projective geometry further augments 
classical projective geometry to restore these concepts. For our purposes, the primary advantage 
of oriented projective geometry is that it allows us to treat an unbounded polyhedron just like 
a bounded polytope — that is, as a convex hull of finitely many vertices (where some of these 
vertices happen to be points at infinity). 

Oriented projective geometry was introduced by Stolfi in the 1980s [37l |38], and has since 
found applications in computer vision |27[ 139] . Kirby presents a leisurely geometric overview in 
[26] , and Boissonnat and Yvinec discuss issues relating to polytopes and polyhedra in |8] . In this 
section we describe the way in which oriented projective geometry extends Euclidean geometry, 
show how we can perform arithmetic in oriented projective geometry using signed homogeneous 
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coordinates, and discuss how unbounded polyhedra in Euclidean space can be reinterpreted as 
projective polytopes. 

Following the original notation of Stolfi, we denote the oriented projective d- space by T^. 
Geometrically, we construct by augmenting the Euclidean d-space as follows: 

• We add a collection of points at infinity. Each ray of passes through one such point, 
and each line of E^ passes through two such points (one in each direction) . Two rays pass 
through the same point at infinity if and only if they are parallel and pointing in the same 
direction. This is illustrated for d = 2 in Figure [TJ 

, . -9 Point at infinity 




"^x) Point at infinity 



Point at infinity 6^ - - " ' 

Figure 1: Euclidean rays and lines passing through points at infinity in T2 

• Travelling beyond these points at infinity, we add a second "hidden" copy of Euclidean 
d-space which we call invisible space; in contrast, we refer to our original copy of E^ 
as visible space. For computational purposes we can ignore invisible space — all of the 
Euclidean structures that we use in this paper will be situated in the visible copy of E^. 
Invisible space is simply a theoretical convenience that endows with a number of useful 
geometric and analytical properties. 



Objects of interest 




Figure 2: Building the oriented projective space T2 

For illustration. Figure [2] shows the full construction of the oriented projective space T2. We 
see from this diagram that T2 has the topological structure of a sphere: the visible points lie on 
the northern hemisphere, the invisible points lie on the southern hemisphere, and the points at 
infinity lie on the equator. This picture generalises to higher dimensions, and Stolfi's original 
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construction shows how points in correspond to points on the unit sphere in M +^ via a 
geometric projection. 

We coordinatise oriented projective c?-space using signed homogeneous coordinates. Each 
point p S is represented by a (d + l)-tuple of numbers, and two {d + l)-tuples x, y S E'^+^ 
represent the same point of if and only if x = Ay for some scalar A > 0. So, for instance, the 
coordinates (2, 3, 1), (4, 6, 2) and (1, |, i) all represent the same point of T2, but the coordinates 
(—2, —3, —1) represent a different point. The coordinates (0, ... ,0) do not represent any point 
of Td at all. 

More specifically, we coordinatise the individual points of as follows: 

• Consider any Euclidean point e = (ei, . . . , ej) G Ed- In T^^, the copy of e in visible space 
is described by the homogeneous coordinates (Aei, . . . , Ae^, A) for any A > 0, and the copy 
of e in invisible space is described by the homogeneous coordinates (— Aei, . . . , — Ae^, — A) 
for any A > 0. 

• Consider any Euclidean ray {/ir | /i > 0} C E^;, where r = (ri, . . . , r^)- In T^, this ray 
in visible space extends to a point at infinity, which is described by the homogeneous 
coordinates (ri, . . . , r^, 0). 

Conversely, any point x £ with homogeneous coordinates (xi, . . . ,Xd+i) can be interpreted 
geometrically as follows: 

• If Xd+i > then x describes the visible Euclidean point (^^, . . . , -^^) G E^. 

• If Xd+i < then x describes the invisible Euclidean point (-^^, . . . , G E^. 

• If Xd+i — then x describes a point at infinity, which is the limiting point of the visible 
Euclidean ray {fj, ■ {xi, . . . , Xd) | A* > 0} C E^;. 

For brevity, we refer to such points as visible, invisible and infinite points of respectively. 
Examples of all three types of point in T2 are illustrated in Figure [3l 

Visible point (-3, 2, 1) 



Infinite point (0, 1, 0) 
Infinite point (—3, —1, 0) — 



Invisible point (—20, —18, 

Figure 3: Different types of points in T2 

A projective hyperplane in is coordinatised hy a {d+ l)-tuple h ~ {hi, . . . , hd+i) G K'^^"'^, 
and contains all points x G for which hixi + . . . + hd+iXd+i = 0, or more concisely h • x = 0. 
Most projective hyperplanes consist of a Euclidean hyperplane in visible space and a Euclidean 
hyperplane in invisible space, joined together by a common set of points at infinity. The only 
exception is the infinite hyperplane described by the {d + l)-tuple (0, . . . , 0, 1), which contains 



-Visible point (4, -2, 1) 



y 



Infinite point (1, —2, 0) 

— Invisible point (6, 1, —1) 
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all of the points at infinity and nothing else. Projective hyperplanes of both types are illustrated 
in Figure m The {d + l)-tuple Ah describes the same projective hyperplane as h for any A ^ 0, 
and the {d + l)-tuple (0, . . . , 0) does not describe any hyperplane at all. 




-Visible point (-2, 3, 1) 



Infinite point (—1, 3, 0) 




-Visible point (1, —6, 1) 





") Projective hyperplane (0, 0, 1) 
(infinite hyperplane) 



Infinite point (-2,-1,0) 



Infinite point (1, — 3, 0) 



Invisible point (3, —6, — 1) 



Projective hyperplane (3, 1, 3) 



Figure 4: Two examples of projective hyperplanes in T2 



A projective half-space in is likewise coordinatised by a (d + l)-tuple h = (/ii, . . . , hd+i) 
G R''+^, and contains all points x G for which h • x > 0. Most projective half-spaces consist 
of a Euclidean half-space in visible space and a Euclidean half-space in invisible space, again 
joined together by a common set of points at infinity. There are two exceptions: the visible 
half-space described by the {d -\- l)-tuple (0, . . . , 0, 1) and the invisible half-space described by 
the {d -\- l)-tuple (0, . . . , 0, —1). The visible half-space contains all visible and infinite points of 
Td, and the invisible half-space contains all invisible and infinite points of T^. All three types 
of projective half-space are illustrated in Figure [5] The {d + l)-tuple Ah describes the same 
projective half-space as h for any A > (but not A < 0), and as usual (0, . . . , 0) does not 
describe any half-space at all. 



For a projective half-space with coordinates h — {hi, . . . , hd+i), the interior of this half- 
space consists of all points x G for which h • x > 0; that is, all points of the half-space that do 
not lie on the boundary hyperplane. It is worth noting that the interior of the visible half-space 
is precisely the set of all visible Euclidean points. 

For any point x G T^; with coordinates (xi, . . . ,Xd+i), the opposite point to x is the point 
with coordinates (— xi, . . . , —x^+i), which we denote ^x. The opposite of a visible point is an 
invisible point, and vice versa. If x is a point at infinity then ^x is another point at infinity, 
and these points lie at opposite ends of the same Euclidean lines. 




An "ordinary" half-space 



The visible half-space 



The invisible half-space 



Figure 5: Different types of projective half-space in T2 
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Arithmetic in is very simple. For distinct and non-opposite points x, y G with coordi- 
nates . . . , Xd+i) and (yi, . . . , y^+i) respectively, the unique line joining x and y consists of 
all points with coordinates (axi + /3yi, . . . , axd+i + Pyd+i) for any a, /? S K, excluding the case 
For the line segment joining x and y we add the constraints a, ^ > 0. In other 
words, the arithmetic of lines in T^; is essentially the same as traditional Euclidean arithmetic in 
R'^+^j except that we discard the constraint a + P — 1 (because signed homogeneous coordinates 
are unaffected by positive scaling). 

We can define convexity in terms of line segments. Consider any set 5' C that lies within 
the interior of some projective half-space. Then S is convex if and only if, for any two points 
x, y S S, the entire line segment joining x and y lies within S. The half-space constraint 
ensures that S cannot contain two opposite points, which means that line segments are always 
well-defined. Note that this half-space constraint is automatically satisfied for any set of visible 
Euclidean points (which all lie in the interior of the visible half-space). 

Consider now some finite set of points x^^^ . . . ,x('^^ S T^;, where each x*^'^ has coordinates 
(xi \ ■ ■ ■ ,x^'^]^), and again suppose that these points all lie within the interior of some pro- 
jective half-space. The convex hull of these points is the smallest convex set containing all 
of x(^\ . . . , x'*"'' . Equivalently, it is the set of all points in with coordinates of the form 
(X] Oiix'f^ , . . . , ^ ctiX^d+i) scalars ai, . . . , > 0, excluding the case ai = . . . = ai = 0. 

It is important to note that, when we restrict our attention to visible Euclidean points, all 
of these concepts — lines, line segments, convexity and convex hulls — have precisely the same 
meaning as in traditional Euclidean geometry. 

As noted earlier, the advantage of oriented projective geometry for us is the way in which 
it handles polytopes and polyhedra. Recall that in Euclidean geometry, a polyhedron (which 
may be bounded or unbounded) can be expressed as an intersection of finitely many half-spaces, 
and that a polytope (that is, a bounded polyhedron) can also be expressed as a convex hull 
of finitely many vertices. Using oriented projective geometry, we are able to express some 
unbounded polyhedra as convex hulls also. The details are as follows. 

We define a projective polytope in to be the convex hull of a finite set of points V ^Td, 
all of which lie within the interior of some projective half-space. If the set V is minimal (that 
is, no point can be removed from V without changing the convex hull), we call the points of V 
the vertices of the projective polytope. 

It is clear that traditional Euclidean polytopes correspond precisely to projective polytopes 
whose vertices are all visible Euclidean points. The situation regarding unbounded Euclidean 
polyhedra is more interesting. 

Suppose that the set ^ C consists of a mixture of visible Euclidean points and points 
at infinity, and (as usual) that all of these points lie within the interior of some projective 
half-space. Then the projective polytope with vertex set V is in fact an unbounded Euclidean 
polyhedron in visible space, "closed off" with additional points at infinity. 

For example. Figure [6] shows a projective polytope in T2, with two visible vertices (2,3, 1) 
and (3,2,1), and two infinite vertices (—1,-1,0) and (0,-1,0). The left-hand diagram shows 
this projective polytope in T2, and the right-hand diagram shows the corresponding unbounded 
polyhedron in E2. 

Conversely, let P be an unbounded Euclidean polyhedron in that lies within the interior 
of some Euclidean half-space. Then P can be expressed as a projective polytope. More precisely, 
there is a projective polytope in that consists of the polyhedron P in visible space "closed 

■^We insist on non-opposite points because there are infinitely many distinct lines in joining x and -ix. 
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Projective polytope Euclidean polyhedron 



Figure 6: A 2-dimensional projective polytope with vertices at infinity 

off" with an extra facet at infinity. The vertices of this projective polytope are the vertices of 
P in visible space with some additional vertices at infinity. 

For a 3-dimensional example, Figure [7] shows a Euclidean polyhedron in E3 with vertices 
(2,0,1), (1,1,2) and (0,2,0), and with facets that extend infinitely far in the negative x, y 
and z directions. The corresponding projective polytope has six vertices: the visible vertices 
(2, 0, 1, 1), (1, 1, 2, 1) and (0, 2, 0, 1), and the three infinite vertices (-1, 0, 0, 0), (0, -1, 0, 0) and 
(0,0,-1,0). 




- (0, 2, 0) 



Figure 7: An unbounded 3-dimensional Euclidean polyhedron 

As in traditional geometry, any projective polytope can be expressed as an intersection of 
finitely many projective half-spaces. Moreover, projective polytopes have the same combinatorial 
properties as Euclidean polytopes, and they support the same combinatorial algorithms (such 
as the double description method for vertex enumeration, which we use later in this paper) H 

•^This equivalence arises because any projective polytope can be transformed into a visible Euclidean polytope 
without changing its combinatorial properties. In the spherical model of Figure[2] the transformation is essentially 
a rotation of the sphere. 
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4 Efficiency-equivalent polyhedra 



Recall that our focus is on generating all efficient extreme outcomes for a multiobjective linear 
program. In other words, our aim is to generate all vertices of the outcome set (a polytope 
in Ep) that are also efficient (that is, not dominated by any other point in Y^). 

One difficulty is that the outcome set Y^ can become extremely complex even for relatively 
small problems [36l 112] , and the efiicient vertices might only constitute a small fraction of the 
set of all vertices of Y^. Several authors address this by working with efficiency-equivalent 
polyhedra — polyhedra in Ep that are typically simpler in structure than Y^ , but whose efficient 
sets precisely match the efficient outcome set Y^. 

There are two such polyhedra that feature prominently in the literature: the bounded poly- 
tope Y^ introduced by Benson [5], and the unbounded polyhedron Y- first described by Dauer 
and Saleh [12l [13] , both of which we describe shortly|3 The goal of this section is to understand 
the structure and complexity of these two polyhedra. 

We begin this section with a general result about efficiency-equivalent polyhedra that justifies 
our focus on vertices of these polyhedra. We then describe the efficiency-equivalent polyhedra 
r° and Y^, and we reinterpret the unbounded polyhedron Y- as a polytope T- using oriented 
projective geometry. The remainder of this section concentrates on the complexity of the Eu- 
clidean polytope Y^ and the projective polytope T-, paying particular attention to the number 
of "unwanted" non-efficient vertices. We find through our analysis that T- is far simpler than 
Y^ (and therefore far superior) in this respect. 

4.1 Definitions and general results 

Our first result shows that any efficiency-equivalent polyhedron for Y^ can be used to generate 
the efficient extreme outcomes for our multiobjective linear program. Ours is a general result; 
Benson [5 and Dauer and Saleh |T2| make similar observations for the specific polyhedra Y^ 
and Y- respectively. 

Lemma 4.1. Let Z C be any efficiency- equivalent polyhedron for Y^ . Then the efficient 
vertices of Z are precisely the efficient vertices of Y^ . That is, the efficient vertices of Z are 
precisely the efficient extreme outcomes for our multiobjective linear program. 

Proof. Let v be an efficient vertex of Y^. By efficiency-equivalence we know that v lies in the 
efficient set Ze- We aim to show that v is a vertex (and thus an efficient vertex) of Z. 

Let F be the smallest-dimensional face of Z containing v. Because Ze is a union of faces 
of Z, every point in F must be efficient. By efficiency-equivalence, the entire face F must lie in 
Y^ also. If V is not a vertex of Z then v lies within the relative interior of the face F, and since 
F C y= it follows that v cannot be a vertex of Y^ , a contradiction. 

Applying the same argument in reverse shows that every efficient vertex of Z is an efficient 
vertex of Y^, and the proof is complete. □ 

We proceed now to define the polyhedra Y^ and Y- , which will occupy our attention for the 
remainder of this section. Both are "extensions" of the outcome set Y^ , in that they enlarge 
the non-efficient portion of Y^ to create a simpler structure. 

*The symbols Y~ and Y— follow Benson's notation; Dauer and Saleh call these Y and Y instead. The symbol 
Y^ is our own notation; Benson simply calls this Y. We avoid the plain symbol Y in this paper because it is 
used with several different meanings throughout the literature. 
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Definition 4.2. We define the Euclidean polyhedron Y~ to be the set of all points in Ep that 
are dominated by any point in Y^. That is, 

Y- = {z G Ep I z < y for some y G Y~ } . 

Choose some point y G Ep that is strictly dominated by all ofY^; in other words, y < y for all 
y G Y^ . We then define the Euclidean polytope Y^ to be 

y'-' = {z G Ep I y < z < y /or some y G } . 




The outcome set Y The efficient outcome set 




The unbounded polyhedron Y- The polytope Y 



Figure 8: Constructing the efficiency-equivalent polyhedra Y- and Y 

Both of these structures are illustrated in Figure |S] (where we use yi, 1/2,1/3 for our three 
coordinate axes to show that we are working in objective space). It is clear that Y- is an 
unbounded polyhedron, and that it can be expressed as the Minkowski sum 

= {y + z |y G r= and z < 0}. (4.1) 

Likewise, it is clear that Y^ is a polytope, since its points are bounded below by y and above 
by Y^ . Geometrically, Y^ truncates the unbounded polyhedron Y~ by intersecting it with a 
large cube that surrounds the original outcome set Y^. 

The following crucial properties of Y^ and Y- are proven by Benson 5 and Daucr and 
Saleh pjj respectively. 

Lemma 4.3. Both Y° and Y^ are efficiency- equivalent polyhedra for Y . In other words, the 
efficient sets Y^, Y^ and Y^ are identical. 
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4.2 The complexity of and 

The unbounded polytope Y- is simpler than in structure, with fewer vertices and a conve- 
nient Minkowski sum representation. However, Y- is more difficult to work with algorithmically 
because it is not a polytope — in particular, it cannot be expressed as the convex hull of its ver- 
tices. This makes it less amenable to traditional algorithms for vertex enumeration and linear 
optimisation. 

We now grant ourselves the best of both worlds by reinterpreting the unbounded polyhedron 
Y- as a convex hull of vertices in oriented projective p-space. 

Definition 4.4. Let T- be the projective polytope in Tp formed as the convex hull of: 
(i) the vertices of the polyhedron Y- in visible Euclidean space; 

(a) the points at infinity defined by the signed homogeneous coordinates (— 1, 0, 0, . . . , 0, 0), 
(0, — 1, 0, . . . , 0, 0), (0, 0, — 1, . . . , 0, 0), (0, 0, 0, . . . , — 1, 0) (there are p such points in 
total). 

In the following series of results, we show that the projective polytope T- is well-defined, 
that it is truly a projective reinterpretation of Y-, and that Definition 14.41 gives the precise 
vertex structure of T~ in oriented projective p-space. 

Lemma 4.5. The projective polytope T- is well-defined. That is, all of the points listed in 
and above lie in the interior of some projective half-space. 

Proof. Because the Euclidean polyhedron Y- has finitely many vertices, there is some a e R 
that satisfies a > vi + . . . -\- Vp for every Euclidean vertex {vi, . . . , Vp) of Y- . When translated 
into oriented projective p-space, it follows that every vertex of Y- lies in the interior of the 
projective half-space with signed homogeneous coordinates (—1, . . . , —1, a). The same is clearly 
true for each infinite point in list ([u]) above, and so the projective half-space (—1, . . . , —1,0) is 
the half-space that we seek. □ 

Lemma 4.6. The visible portion of the projective polytope T- is precisely the Euclidean poly- 
hedron Y- . 

Proof. We prove this by first showing that every visible Euclidean point of T- belongs to 
and then that every point of the visible Euclidean polyhedron Y- belongs to T-. 

• Let t G Tp be some visible Euclidean point in T-. We first observe that any convex 
combination of points in list ([!]) above is a visible Euclidean point in 1"-, and that any 
convex combination of points in list ([n]) above is a point at infinity with signed homoge- 
neous coordinates (— . . . , — /3p, 0) for some scalars /3i, . . . , /3p > 0. Since t is a convex 
combination of points in both lists, it follows that t has signed homogeneous coordinates 
of the form {^yi — \ . . . , iiyp — \f3p, where /i. A, /3i, . . . , /3p > and where (yi, . . . , yp) 
are the Euclidean coordinates of some point uiY- . 

Because t is visible we must have /i > 0, so we can rescale the coordinates of t to 
{yi — A/3i//i, . . . ,yp — Xf5p/ II, 1). It is clear then that t is a visible point whose Euclidean 
coordinates are formed by adding the non-positive vector — A//i(/3i, . . . , f3p) to some point 
in Y-. From equation (|4.1I) it follows that t is a visible Euclidean point in Y~. 
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• Now let y £ Tp be a point of in visible Euclidean space with signed homogeneous 
coordinates (yi, . . . , 1). Let z = (zi, . . . , Zp, 1) be an efficient point of Y- that domi- 
nates y. We can express the coordinates of y as {zi — Pi, . . . , Zp — fHp, 1) for some scalars 
Pit ■ ■ tPp > 0, which shows that y is a convex combination of the efRcient point z with 
the points at infinity in list (|ii]) above. 

Recalling that the efficient set is a union of faces of Y-, let F be some efficient face 
of Y- containing z. Because Y- is efficiency-equivalent to the bounded polytope F", 
the efficient face F must likewise be bounded. Therefore z can be expressed as a convex 
combination of vertices of F, which are of course vertices of Y- also. This allows us to 
express y as a convex combination of points in lists Q and and so y belongs to the 
projective polytope T-. □ 

Lemma 4.7. The vertices ofT- are precisely the points listed in Definition \4-4\ In other words, 
no point can be removed from these lists without changing the convex hull. 

Proof. To prove this we must show that no point in lists Q and (|n]) from Definition 14.41 can 
be expressed as a convex combination of the others. Suppose that list (P contains the points 
Vi, . . . , Vt (these are the vertices of Y-), and that list (jn]) contains the points fi, . . . , fp (these 
are all points at infinity). 

• No ii can be expressed as a convex combination of the other points in these lists, since 
each other point introduces some non-zero coordinate that is zero in f^. 

• Suppose that some can be expressed as the convex combination J2j^i ^j^i + X!fc7fcffc 
where each 0^,7^ > (as usual we work in homogeneous coordinates). Some Uj must be 
non-zero, since otherwise v.; would be a point at infinity. Some 7^ must also be non-zero, 
since we cannot express a vertex of Y- as a convex combination of other vertices. 

It follows that we can express in the form (yi— . . . , ?/p— /3p, 1), where y = {yi, . . . , j/p, 1) 
is a convex combination of vertices oiY- (in particular, y G Y~), and where /3i , . . . , /3p > 
with at least one (3k strictly positive. This places in the interior of the line segment 
joining y with the point (yi — 2/3i, . . . ,yp — 2/3p, 1). However, both endpoints of this line 
segment lie in the polyhedron Y-, which is impossible since is a vertex of Y-. □ 

We can draw some important conclusions from these results. The first is that T- is precisely 
the efficiency-equivalent polyhedron Y-, recast as a polytope using oriented projective geometry. 
Furthermore, Lemma 14.71 gives an insight into the complexity of this polytope, as seen in the 
following result. 

Corollary 4.8. The visible Euclidean vertices ofT- are precisely the efficient extreme outcomes 
for our multiobjective linear program. That is, the only vertices of T- that are not efficient 
extreme outcomes are the p vertices at infinity. 

Proof. From Lemma 14. 7[ the visible Euclidean vertices of T- are precisely the vertices of the 
unbounded polyhedron Y- C Ep. Moreover, Lemma |4T] shows that the efficient vertices of Y- 
are precisely the efficient extreme outcomes. 

Suppose that Y- has some non-efficient vertex v £ Y-. Then v must be dominated by 
some other point x £ Y-, so that the difference d = x — v is non-negative and non-zero. Now 
construct the point x' = v — d. It is clear that x' £ Y- and that v lies strictly inside the line 
segment [x, x'] C Y-, contradicting the claim that v is a vertex of Y-. □ 
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It follows from Corollary 14.81 that the projective polytope T- gives an extremely compact 
representation of the efficient extreme outcomes — the number of additional "unwanted" vertices 
grows linearly with the number of objectives p, and does not depend at all on the problem size 
(n, m). 



We turn now to study the complexity of the Euclidean polytope Y . We begin with a complete 
categorisation of vertices of F^, and we follow with a discussion of what this means for best- 
case and worst-case complexity. Throughout the remainder of this section we work entirely in 
Euclidean coordinates. 

Definition 4.9. For each set of indices S C {1, . . . we define the set of points V{S) C Ep 
as follows. Denote the p objectives by yi, . . . , yp, and consider the smaller multiobjective linear 
program obtained by considering only the objectives {yi \ i G S} (in other words, by deleting the 
ith row of the objective matrix for each i ^ S). This smaller problem yields a smaller set of 
efficient extreme outcomes in the space ]E|s|- 

We extend these efficient extreme outcomes to points in Ep as follows. For each deleted 
objective i ^ S, we fix the coordinate yi to the value % (using the point y G Ep that we describe 
in Definition \4-.^ . The resulting set of points in Ep is denoted V{S). 

If S is the empty set then there is no linear program to solve, and we simply define V(Jl)) — 



It is worth noting that for the full set S — {1, . . . ,p}, the points in V{S) are simply the 
efRcient extreme outcomes for our original (full-dimensional) multiobjective linear program. 

The various sets V{S) are illustrated in Figure [9l which continues the earlier example of 
Figure m We see in this diagram that the sets V{S) together describe all of the vertices of Y^. 
This is no accident, as shown by the following result. 

Theorem 4.10. The sets V{S) are pairwise disjoint; that is, if S ^ S' then V{S)r\V{S') = 0. 
Moreover, the union of all such sets V{S) over all S C {1, . . . ,p} gives the complete vertex set 
for the polytope Y'^ . 

Proof. Pairwise disjointness is simple to prove. Because y is strictly dominated by all y e Y^ , 
we find that for any v = {vi, . . . ,Vp) e V{S) we have Vi — yi if and only if z e S*. It is then clear 
that V cannot belong to both V{S) and V{S') for distinct sets S ^ S' . 



4.3 The complexity of Y 



{?}■ 




Figure 9: The sets V{S) in relation to the polytope Y 
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It remains to show that every eleraent of each set V{S) is a vertex of Y , and that every 
vertex of can be expressed as an element of some V{S). For this it is convenient to define the 
hnear projection Trg: Ep E|s| that simply deletes all coordinates yi for i ^ S. For example, 
if S' = {1, . . . ,p} then ns is the identity map, and if S" = {i} then its merely extracts the «th 
coordinate. 

From Definition 14.91 we have v g V{S) if and only if (i) Vi — yi for each i ^ S, and (ii) the 
projection 7r5(v) is an efhcient vertex of the polytope tts(Y^). However, applying Lemma 14.31 
to the smaller multiobjective linear program with objectives {yi \ i G S} shows that nsiY^) is 
efficiency- equivalent to Trs{Y^). This means that v G V{S) if and only if (i) Vi — yi for each 
i ^ S, and (ii') 7r5(v) is an efficient vertex of irsiY^). 

We can now conclude the proof as follows: 

• Suppose that v e V{S) for some set S* ^ 0. We aim to show that v is a vertex of Y^. If 
not, we can express v = i(u + w) for some distinct u,w G Y^. We proceed to derive a 
contradiction by showing that v = u = w. 

— For any i ^ S we have u,; = yi by definition of V{S), and Ui,Wi > yi by definition of 
y^. Since v = i(u + w) it follows that Vi = Ui = Wi. 

— By linearity, 7r5(v) = i (7r5(u) -l-7r5(w)) where 7r5(u), 7r5(w) G iisiY^). However, 
because v G V{S) we know that 7rs(v) is a vertex of irs^Y^), which is only possible 
if 7rs(v) = 7rs(u) = 7r5(w). That is, Vi = Ui = Wi for all z G 5. 

This establishes v = u = w, and we conclude that v must indeed be a vertex of Y^. 

We handle the case 5 = by observing that V{%) contains the single point y, and that y 
is a vertex of Y^ because it uniquely minimises the linear functional y >— > j/i + • . ■ + 2/p. 

• Now let V be some vertex of Y^, and let S = {i\vi ^ yi}. We aim to show that v G V{S). 
If 5 = then v y G 1^(0), so assume that S is non-empty. We already have Vi ~ % for 
each i ^ S, so we only need to show that 7r5(v) is an efficient vertex of tts{Y^). 

— If TTsiy) is not a vertex of tts{Y^), there must be some u,w G Y^ for which 
T^'si'v) = f ('''■s(u) -I- 7rs(w)) but 7rs(u) ^ 7r5(w). Construct the points u',w' from 
u, w respectively by replacing each ith coordinate with yi for all i S. 

It is clear that u', w' G Y^ (since y <\i' <\i and y < w' < w). Moreover, for i ^ 5 
we have Vi — u'i = w'i, and for i G S we have Vi = ^{u'i + w'^). Finally, we note that 
7rs(u) 7^ 7rs(w) implies that u' ^ w'. We therefore have v = i(u' + w') but u' w', 
contradicting the assumption that v is a vertex of Y^ . 

— If 7rs(v) is not efficient in irsiY^), there must be some u G Y^ for which 7r5(u) > 
7''s(v) but 7rs(u) 7rs{v). Again, construct u' from u by replacing the ith coordinate 
with yi for each i ^ S. This time we have u' G Y^, u' ^ v, and the relations 
yi = Vi ~ u'i for i ^ S and yi < Vi < for i <E S. 

Let z = V — e(u' — v) for some small e > that ensures y < z (by the relations above, 
such an e exists). Then z G Y^ and v is in the interior of the line segment [z, u'], 
again contradicting the assumption that v is a vertex of Y^ . 

Together these observations show that 7rs(v) is an efficient vertex of nsiY^), whereupon 
the proof is complete. □ 
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Theorem 14 . 1 01 has important imphcations for the complexity of the polytope Y^. Recall that 
all of the efficient vertices of are held in the single set • ■ ■ ,p})- This means that for 

every proper subset S* C {1, . . . the corresponding V(S) is filled with non-efficient vertices 
ofr°. This can yield an extremely large number of non-efficient vertices — each V{S) describes 
the efficient extreme outcomes for some other multiobjective linear program with up to p — 1 
objectives, and we already know that the number of such outcomes can grow very large even 
for relatively small values of p, n and m [321 112] • 

This argument shows intuitively that Y^ can have an extremely large number of non-efficient 
vertices. We follow now with two corollaries of Theorem 14.101 that back this up with explicit 
formulae. The first result shows that, even in the best case scenario where Y^ is as simple as 
possible, we still obtain an exponentially large number of non-efficient vertices. 

Corollary 4.11. The polytope Y^ always has at least 2^ — 1 non-efficient vertices; that is, 2^ — 1 
vertices that are not efficient extreme outcomes for our original multiobjective linear program. 

Proof. As noted above, every proper subset S C {1, . . . ,p} yields a corresponding V{S) that 
consists entirely of non-efficient vertices. There are 2^* — 1 proper subsets of {1, . . . and so 
we have at least 2^ — 1 non-cfficicnt vertices in total. □ 

The corollary above shows that the growth rate of non-efficient vertices relative to the number 
of objectives p is unavoidably exponential. In our second corollary, we show that the number of 
non-efficient vertices can also grow at a severe rate relative to the problem size [n,m). This is 
an asymptotic result, and we use the standard notation for complexity whereby O(-) and ri(-) 
denote asymptotic upper bounds and lower bounds respectively. 

Corollary 4.12. Suppose we fix p and allow the problem size {n,m) to vary. Then there are 
multiobjective linear programs for which the number of non-efficient vertices of Y^ grows at 
least as fast as f7(nL(P-i)/2J ). 

Proof. We establish this result using dual cyclic polytopes. For any integers 2 < d < k, the (d, k) 
dual cyclic polytope is a d-dimensional polytope with precisely k facets and 

vertices. For further information on dual cyclic polytopes the reader is referred to a standard 
reference such as Griinbaum |22j. 

For us, the important fact about dual cyclic polytopes is that for a fixed dimension d, 
the number of vertices grows at least as fast as rj(fcL<i/2J) This is because the sum of binomial 
coefficients in equation (|4.2p can be rewritten as a [(i/2j -degree polynomial in k (taking separate 
cases for odd and even d). 

Our strategy then is as follows. Suppose we are given some number of objectives p and some 
number of variables n. We assume here that 3 < p < n (since we are interested in the growth 
rate as n becomes large, and for p < 2 the result is trivial). Let H be the hyperplane in Ep 
defined by the equation yi + ■ . . + yp = 1. We now embed a {p — ^,n) dual cyclic polytope P 
within the hyperplane H (a straightforward procedure since dimH — p ~ 1). 

The polytope P C Ep can be expressed as the outcome set Y^ for some multiobjective linear 
program with p objectives and n variables (its n facets arise from the inequalities x > in the 
solution space E„). Although P is not full-dimensional in Ep, this is just a convenience to keep 
our arithmetic simple. If a full-dimensional example is ever required then we simply extend P 
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to a p-dimensional cone; the analysis becomes messier but the final growth rate of fi(nL(P-i)/2J ) 
remains the same. 

We can now make the following observation. 

Claim: Every vertex ofY^ (that is, P) gives an efficient extreme outcome for some 
smaller multiobjective linear program that uses only p — 1 of our p objectives. 

We will prove this claim shortly. In the meantime, we examine its consequences. By Theo- 
rem IHTOl it follows that each vertex of P corresponds to some element of V{S) where S is some 
{p — l)-element subset of {1, . . . ,p}. Moreover, because P lies entirely within the hyperplane 
^ — 1, any two vertices of P must differ in at least two coordinates. It follows that any two 
distinct vertices of P must correspond to either two distinct elements of the same set V{S), or 
else elements of two distinct sets V{S), V{S'). 

As a result, the sum of jl^CS*)! over all (p — l)-element sets S is at least the number of vertices 
of P. However, each such S* is a proper subset of {1, ... and we can therefore deduce that 
the total number of non-efficient vertices of is at least the number of vertices of P. The 
number of vertices of P is described by equation (|4.2|) . which we noted earlier grows at least as 
fast as fi(fcL'^/2J) = fi(nL(P-i)/2J ). Thus our final result is established. 

All that remains then is to prove the claim above. Let v be any vertex of Y^. It follows 
that V is the unique intersection of Y^ with some supporting hyperplane C Ep. Suppose 
that H is defined by the equation aiyi + . . . + ctpj/p = 0, where Y^ lies in the closed half-space 
aiyi -I- ... -I- apyp < 0. 

We can "tilt" this hyperplane as follows. Let amin = min{ai}, and consider instead the 
hyperplane H' C Ep defined by the equation (ai — amin)yi + . . . + (ctp — amin)yp = 0. Because 
every point in Y^ satisfies '^yi = 1, we find that H' is also a supporting hyperplane for Y^, 
that Y^ lies in the closed half-space (ai — amin)2/i -I- ... -I- (ctp — aynin)yp < 0, and that v is the 
unique intersection of Y^ with H' . 

More importantly, at least one of the coefficients {ai — amin) is zero, and all of the remaining 
coefficients of H' are non-negative. It follows that (ai — Q;inin)yi -I- ... -I- (ctp — cuminjyp is a 
linear functional with non-negative coefficients in at most p — I objectives, and that v is the 
unique maximum for this linear functional in the outcome set Y^ . As a result, v indeed gives 
an efficient extreme outcome for some smaller multiobjective linear program using only p — 1 oi 
our p objectives. □ 

5 The outer approximation algorithm 

In this section we illustrate the power of our techniques by applying them to the outer approxi- 
mation algorithm, which enumerates efficient extreme outcomes for a multiobjective linear pro- 
gram. Introduced by Benson |5j, this algorithm iteratively constructs the efficiency-equivalent 
polytope Y^ using methods from linear programming, univariate search techniques and poly- 
tope vertex enumeration algorithms. The algorithm is designed to work with polytopes, not 
unbounded polyhedra, which explains the choice of Y^ as the final target. 

We begin in Section 15.11 by describing the original outer approximation algorithm in detail. 
In Section r5.2l we develop a new outer approximation algorithm by working in oriented projective 
geometry, which allows us to replace the target polytope Y^ with the much simpler projective 
polytope T-. In Section [5.31 we examine how these changes affect the runtime complexity. 
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5.1 The original algorithm 



The algorithm originally described by Benson begins with a simplex in Ep that completely sur- 
rounds the outcome set F= , and then successively truncates this simplex along hyperplanes until 
it is eventually whittled down to the target polytope Y^. This overall procedure is illustrated 
in Figure [TUl (which continues the earlier example from Figure [5]). The details are as follows. 




Figure 10: Stepping through the original outer approximation algorithm 

Algorithm 5.1 (Outer approximation algorithm). This algorithm builds by constructing 
a series of polytopes Pq ^ Pi ^ ■ ■ ■ ^ Pk = Y^ in Ep (where the number of polytopes k is 
not known in advance). Each polytope is stored using two representations: as a convex hull of 
vertices in Ep, and as an intersection of half- spaces in Ep. 

To initialise the algorithm, choose an arbitrary point y in the interior of Y^ , and set the 
initial polytope Pq to be a simplex with vertex set {y, y + aei, . . . , y + aep}. Here y is the 
point chosen in Definition \4-.S\ each represents the ith unit vector in Ep, and the constant a 
is chosen so that the simplex contains the entire outcome set Y^ . 

Now iteratively construct each polytope Pi+i from the previous polytope Pi as follows: 

Step 1. Search through the vertices of Pi for one that does not belong to the target polytope Y^ . 
If no such vertex exists (which indicates that Pi = Y^), then output the efficient vertices 
of Pi and terminate the algorithm. 

Step 2. Otherwise let v be such a vertex, and compute the unique point x G Ep at which the line 
segment joining v with y crosses the boundary of the target polytope Y^ . 

Step 3. Compute a half-space H C Kp that contains the target polytope Y'^ , and that has the 
point X in its bounding hyperplane. 

Step 4. The next polytope is now Pi+i = Pi D H . Compute a half-space representation of Pi+i 
(which is trivial) and a convex hull representation (which is not). 
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There are significant details to be filled in for each of these steps, for which we refer the 
reader to Benson's original paper [S]. In brief: 

• In the initialisation phase, both the constant a and the interior point y G can be 
found by solving a linear program, and the half-space representation of the initial simplex 
Pq is then simple to compute. 

• In Step [1] of the iteration, we can test for membership in by testing the feasibility of 
another linear program. The boundary point x in Step [2] is found by combining feasibility 
tests with standard univariate search techniques (Benson uses a simple bisection search), 
and the half-space H in Step |3] uses the solution to a dual linear program. 

• In the final output phase, the efficient vertices of Pi — Y^ are precisely those vertices that 
are strictly dominated by y. 

The only step not covered in the list above is Step 01 which requires us to compute half- 
space and convex hull representations for the new polytope Pi+i- The half-space representation 
is simple: just append H to the list of half-spaces for the previous polytope Pi. The convex hull 
representation is more complex, and requires iterative vertex enumeration techniques. These 
are typically based on the double description method of Motzkin et al. [29j , though the central 
idea has been rediscovered by many different authors and goes under a variety of names. See 
Chen et al. [S] and Horst et al. [53] for alternative treatments. 

The fundamental principle behind the double description method is the following result: 

Theorem 5.2. Let P QEp be a polytope, and let H QKp be a half-space. Then the vertices of 
the polytope P C] H can be obtained from the vertices of P using the following procedure. 

First partition the vertices of P into three sets S^, Sq and according to whether each 
vertex lies within the interior of H , on the bounding hyperplane of H , or external to H respec- 
tively. 

Each vertex in Sq or S+ then becomes a vertex of P H H . Furthermore, for each pair of 
vertices u £ S+ and v € that are adjacent in P, the point where the line segment [u, v] cuts 
the bounding hyperplane of H also becomes a vertex of PClH. The polytope PHH has no other 
vertices besides those described here. 



H 




Figure 11: Intersecting the polytope P with the half-space H 

This procedure is illustrated in Figure [TT] The only complication of Theorem 15.21 is de- 
termining which pairs of vertices u G 5+ and v G S*- are adjacent (i.e., joined by an edge) 
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in the original polytope P. There are many different adjacency tests in the Uterature; in this 
paper we use the combinatorial adjacency test, which is found to be fast and simple in practice 
[T5] . and which translates well to the setting of oriented projective geometry (as seen later in 
Theorem [5H). 

Lemma 5.3. Let P C Ep be a polytope with half-space representation P = Hi H H2 H . . . n Hk, 
where each Hk is a half-space in Ep. Then two vertices u, v of P are adjacent if and only if 
there is no other vertex z of P with the following property: 

• For each half-space Hi, if both u and v lie on the bounding hyperplane of Hi then z lies 
on the bounding hyperplane also. 

For proofs of Theorem l5.2l and Lemma l5.3l as well as further details on the double description 
method, the reader is referred to Fukuda and Prodon [18) . 

5.2 A new algorithm 

We now apply the results of Section to the outer approximation algorithm, yielding a new 
algorithm with significant benefits over the original. Our strategy is to work in the oriented 
projective space Tp instead of the Euclidean space Ep. This allows us to iteratively construct the 
projective polytope T- instead of the Euclidean polytope Y^, thereby removing an exponential 
number of non-efficient vertices that would otherwise clutter up our computations. We discuss 
the running time benefits in more detail in Section 15.31 

As with the original algorithm, the final output for our new algorithm is the set of all efficient 
extreme outcomes for our multiobjective linear program. Because we work in oriented projective 
space, we perform all calculations using signed homogeneous coordinates. The details are as 
follows. 

Algorithm 5.4 (New outer approximation algorithm). In this algorithm we build the projective 
polytope T- by constructing a series of projective polytopes Qo ^ Qi ^ ■ ■ ■ ^ Qk — T- in Tp. 
As usual each polytope is stored using two representations: as a convex hull of vertices in Tp, 
and as an intersection of projective half-spaces in Tp. 

To initialise the algorithm, construct the point y™^'' with signed homogeneous coordinates 
iuT^^^ ■ ■ ■ ^Vp^^^T^)' where each y™3,x individually maximises the ith objective. Choose an ar- 
bitrary point y in the interior of T-, and set the initial polytope Qo to the simplex with one 
visible Euclidean vertex y™'" and p vertices at infinity with signed homogeneous coordinates 
(-1,0,0,. ..,0,0), (0,-1,0,. ..,0,0), (0,0,0,. ..,-1,0). 

Now iteratively construct each polytope Qi+i from the previous polytope Qi as follows: 

Step 1. Search through the visible Euclidean vertices of Qi for one that does not belong to the 
target polytope T- . If no such vertex exists (indicating that Qi = T-), then output all 
visible Euclidean vertices of Qi and terminate the algorithm. 

Step 2. Otherwise let v be such a vertex, and compute the unigue point x GTp at which the line 
segment joining v with y crosses the boundary of the target polytope T- . 

Step 3. Compute a projective half-space H C Tp that contains the target polytope T- , and that 
has the point x in its bounding hyperplane. 

Step 4- The next polytope is now Qi+i — Qi Cl H . Compute a half-space representation and a 
convex hull representation of Qi+i . 
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Figure 12: The new algorithm as seen from visible EucHdean space 



The overall procedure is illustrated in Figure [T^ (which only shows the portions of Qo, Qij • ■ • 
in visible Euclidean 3-space, and does not show the additional points at infinity). It should be 
noted that all of the geometric entities used in this algorithm lie in the interior of a projective 
half-spaceH which means that all of the relevant concepts (such as polytopes, convex hulls and 
line segments) are properly defined. 

We will shortly discuss the implementation of the various steps of this algorithm. However, 
our most immediate task is to prove the correctness of the algorithm as a whole. 

Theorem 5.5. Algorithm \5.4\ is correct (i.e., it outputs precisely the efficient extreme outcomes 
for our multiobjective linear program), and it terminates in finitely many iterations. 

Proof. We prove the correctness of this algorithm in several stages. 

(i) For each polytope Qi we have Qi D T- ; 

This is simple to prove by induction. The visible Euclidean points of the initial simplex Qo 
are precisely those with signed homogeneous coordinates of the form {y^'^^ — f3i, . . . , y™^^ — 
Pp, 1) for Pi, . . . , (3p > 0, which includes all visible Euclidean points of T-. Moreover, the 
vertices at infinity for Qq are precisely the vertices at infinity for T- fLemma 14. 7p . Thus 
Qo 3 T^. 

For the inductive step, we observe that Qi+i is the intersection of Qi with a half-space 
H (Step[3l). It follows that if Qi D then Q,+i D also. 

(ii) For each polytope Qi, the vertices at infinity are precisely the points (— 1, 0, 0, . . . , 0, 0), 
(0,-1,0,. ..,0,0), (0,0,0,..., -1,0).- 

Let Hoo denote the infinite hyperplane in Tp (as defined earlier in Section [3]). The in- 
tersection of the initial simplex Qq with H^o is precisely the (p — l)-simplex whose ver- 
tices are (-1, 0, 0, . . . , 0, 0), (0, -1, 0, . . . , 0, 0), (0, 0, 0, . . . , -1, 0). However, Defini- 
tion 23] shows that the intersection of T- with H^o is this same {p — l)-simplex. Because 
Qo D Qi 3 . . . 3 T-, it follows that Qi D H^o = T- n H^o for each i, and in particular the 
vertices of Qi at infinity are the same p points described above. 

^One such half-space is {(xi, . . . , Xp+i) G Tp \ xx + . . . + Xp < {yf''^ + . . . + y^""^ + 1) ■ Xp+i } . 
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(iii) The termination condition is correct, i.e., we terminate if and only if Qi = 

Step [T] instructs us to terminate if and only if no visible Euclidean vertex of Qi lies outside 
T- . By fact ([u]) above, this is equivalent to the condition that no vertex of Qi lies outside 
T-. Because Qi and T- are both projective polytopes with Qi D T-, this is equivalent 
to the condition that Qi = T-. 

(iv) The output is correct, i.e., we output precisely the set of efficient extreme outcomes for our 
multiobjective linear program: 

By fact (pn)) above, if we terminate then we have Qi = T-, whereupon Step [T] instructs us 
to output all visible Euclidean vertices of T-. Corollarv 14.81 shows that these are precisely 
the efficient extreme outcomes for our multiobjective linear program. 

Together these facts show that, if the algorithm terminates, it produces the correct results. 

Our final task is to prove that the algorithm terminates in finitely many steps. Let dQi denote 
the boundary of the projective polytope Qi. Because Qi 3 T-, the intersection dQi DT- must 
be a union of faces of T- . 

We observe that each intersection dQi+iCiT- is strictly larger than dQiCiT-, since it contains 
the new point x identified in StepO In particular, we know that x ^ dQi, since x lies strictly 
between the vertex v and the interior point y. However, we also know that x € dQi+i n T-, 
since Step O ensures that x G T- and Step [3] ensures that x G dQi+i. 

It follows that each intersection dQi+i H T- contains strictly more faces of T- than the 
previous intersection dQi n T-. Because T- has finitely many faces, this sequence cannot grow 
forever and the algorithm must terminate. □ 

The individual steps of Algorithm l5.4l are performed in much the same way as for the original 
algorithm. 

• In the initialisation phase, the point y™'^'^ is straightforward to compute, and for the inte- 
rior point y G T- we can simply use y from Definition 14.21 The half-space representation 
of the initial simplex Qo involves the p "ordinary" half-spaces with signed homogeneous 
coordinates (-1, 0, 0, . . . , 0, 2/?-"), (0, -1,0,..., 0, yf^), . . . , (0, 0, 0, . . . , -1, y'^'^^), as wdl 
as the visible half-space (0, 0, 0, . . . , 0, 1). 

• In Step [1] we are required to test whether a visible Euclidean vertex belongs to For this 
we can restrict our attention to visible Euclidean space (where T- becomes the Euclidean 
polyhedron Y-), and we can employ precisely the same feasibility test used in the original 
algorithm. 

• Likewise, Steps [2] and [3] can be performed in visible Euclidean space using the same feasi- 
bility tests, univariate search techniques and dual linear programs as before. In Step [3] it 
is trivial to convert the resulting half-space H back into oriented projective p-space. 

This leaves Step HI where we must enumerate the vertices of the new polytope Qi+i. For 
this we can use the double description method directly in oriented projective p-space: 

Theorem 5.6. The results of Theorem \5.S\ and Lemma \5.3\ also hold true in oriented projective 
p- space. That is, they remain true if we replace all occurrences ofEp, polytopes, half-spaces 
and hyperplanes with Tp, projective polytopes, projective half-spaces and projective hyperplanes 
respectively. 
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Proof. Instead of offering a direct proof, we can transform these results into their Euchdean 
counterparts using a powerful tool known as a projective mapping. Projective mappings are 
geometric transformations in Tp that preserve collinearity and concurrency. They map projective 
hyperplanes to projective hyperplanes and projective half-spaces to projective half-spaces (and 
therefore projective polytopes to projective polytopes). Every projective mapping has an inverse 
projective mapping, and for any projective half-space H QTp there is a projective mapping that 
takes H to the visible half-space. For more information on projective mappings, the reader is 
referred to Boissonnat and Yvinec [5]. 

Using projective mappings, the proof of Theorem 15.61 becomes simple. If P C Tp is a 
projective polytope then it lies in the interior of some projective half-space iJ C Tp. Let tt be 
a projective mapping that maps H to the visible half-space. The interior points of H map to 
visible Euclidean points, and P therefore maps to a visible Euclidean polytope. 

We can now apply Theorem 15.21 and Lemma [5.31 directlv to the Euclidean polytope 7r(P). 
Because both results are expressed purely in terms of lines, intersections, half-spaces and 
hyperplanes — all notions that are preserved under a projective mapping — the corresponding 
results are therefore true for the original projective polytope P. □ 

When using the double description method in this way to enumerate the vertices of Qi+i, it 
is important to take into account all of the vertices of Qi, including the vertices at infinity. 

This final point highlights a key motivation for working in oriented projective geometry. 
In Euclidean space the double description method works only for polytopes, not unbounded 
polyhedra, because it requires a description of the polytope as a convex hull of its vertices. 
Oriented projective geometry allows us to use the double description method with unbounded 
Euclidean polyhedra (such as Y-) by treating them as polytopes — that is, convex hulls of 
vertices — in Tp. 

5.3 Running time 

The primary benefit of our new outer approximation algorithm is that we avoid constructing 
and analysing a large number of "unnecessary" vertices that are not efficient extreme outcomes 
for our multiobjective linear program. As shown in Section |4l the projective polytope T- (used 
in the new algorithm) has only p unnecessary vertices, which are the p vertices at infinity. In 
contrast, the Euclidean polytope (used in the original algorithm) has a significant number of 
unnecessary vertices — their number is always exponential in the number of objectives p (Corol- 
larv 14.111) . and can also grow at a severe rate relative to the problem size (n, m) (Corollarv l4.12p . 

Just as has significantly more vertices than T- , we also expect the intermediate polytopes 
Pi, P2, . . . from the original algorithm to have significantly more vertices than Qi, Q2, . . . from 
the new algorithm (recall that these families of polytopes successively approximate Y^ and T- 
respectively) . This growth in the number of vertices has a direct impact on the running time of 
the outer approximation algorithm (Algorithm 15. ip : 

• In Step [2 the time required to search for a vertex v outside the target polytope Y^ is 
directly proportional to the number of vertices in the intermediate polytope Pi. 

• In Step m the time required to run the double description method is proportional to the 
cube of the number of vertices: if the intermediate polytope Pj has V vertices, we have 
0{V^) potential pairs u G 5*+ and v G 5_ (Theorem 15. 2p . and for each such pair the 
adjacency test requires us to search through all of the vertices again (Lemma l5.3p . 
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• More generally, the number of stages in the algorithm as a whole (that is, the number of 
intermediate polytopes Pi, . . . ,Pfc) depends on the number of vertices: each stage corre- 
sponds to some intermediate vertex v outside the target polytope Y^, and so with more 
intermediate vertices we expect more intermediate polytopes in total. 

All of these observations indicate that the new algorithm should offer significant benefits 
over the original algorithm in terms of running time. 

The only added complication for the new algorithm is that we must perform arithmetic in 
oriented projective p-space, where we work with p + 1 coordinates at a time instead of p. The 
cost of this "extra coordinate" is negligible, and (unlike the number of vertices) does not affect 
the asymptotic growth of the running time at all. 

6 Conclusions 

In this paper we show how oriented projective geometry can be applied to the field of multiob- 
jective linear programming, allowing us to work with efficiency-equivalent projective polytopes 
that are significantly less complex than their Euclidean counterparts. Polytopes (as opposed 
to unbounded polyhedra) are important because they support techniques such as the double 
description method, which operates on convex hulls of finitely many vertices. The reduction in 
complexity is important because it leads to a significant reduction in running time. 

As a concrete illustration, we apply these techniques to the outer approximation algorithm, 
which is used to generate all efficient extreme outcomes for a multiobjective linear program. By 
working in oriented projective geometry we obtain a considerable reduction in the number of 
non-efficient vertices that the algorithm generates, and we show how this directly benefits the 
running time as a result. 

We support these results with explicit asymptotic estimates on the number of non-efficient 
vertices generated in the old Euclidean setting and the new oriented projective setting. Using 
oriented projective geometry, the number of non-efficient vertices is precisely the number of 
objectives p. In the old Euclidean setting this number is at least 2^ — 1, and it can grow as fast 
asl7(nL(p-i)/2J). 

For applications with even a moderate number of objectives, this growth rate can be ex- 
tremely high. For example, Steuer describes a blending problem with p = 8 objectives [34 ; here 
the number of non-efficient vertices in the old algorithm can grow like il{n^), which for large 
problem sizes is a significant burden. When the number of objectives is higher, as in the aircraft 
control design problem of Schy and Giesy with p — 70 objectives [35], even the best-case growth 
rate of 2^ — I becomes crippling. In contrast, the new techniques in this paper deliver just 8 
and 70 non-efficient vertices respectively, which is negligible in comparison. 

Although we focus on the outer approximation algorithm for our application in Section [5l 
the underlying techniques of Sections [3] and |4] are quite general. We expect that the methods 
described in this paper can find many fruitful applications throughout the field of multiobjec- 
tive linear programming, and that with further research the geometric and analytical insights 
described here can lead to new discoveries in other areas of high-dimensional optimisation. 
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